Graphene - a nearly perfect fluid 
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Hydrodynamics and collision dominated transport are crucial to understand the slow dynamics of 
many correlated quantum liquids. The ratio r]/s of the shear viscosity n to the entropy density s is 
uniquely suited to determine how strongly the excitations in a quantum fluid interact. We determine 
77/s in clean undoped graphene using a quantum kinetic theory. As a result of the quantum criticality 
of this system the ratio is smaller than in many other correlated quantum liquids and, interestingly, 
comes close to a lower bound conjectured in the context of the quark gluon plasma. We discuss 
possible consequences of the low viscosity, including pre-turbulent current flow. 
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Graphene [TJ |2] , attracts a lot of attention due to the 
massless relativistic dispersion of its quasiparticles and 
their high mobility. Recently, it was shown that this 
material offers a unique opportunity to observe trans- 
port properties of a plasma of ultrarelativistic particles 
at moderately high temperatures [3 . Undoped graphene 
is located at a special point in parameter space where 
the Fermi surface shrinks to two points, and in many re- 
spects it behaves similarly as other systems close to more 
complex quantum critical points 0]. Due to its massless 
Dirac particles graphene also shares interesting proper- 
ties with the ultrarelativistic quark gluon plasma. The 
latter, surprisingly, has an unexpectedly low shear vis- 
cosity, as was observed in the dense matter balls created 
at the relativistic heavy ion collider RHIC [5]. We show 
here that an analogous property can be found in undoped 
graphene, reflecting its quantum criticality. 

The shear viscosity 77 measures the resistance of a fluid 
to establishing transverse velocity gradients, see Fig. [T| 
The smaller the viscosity, the higher the tendency to tur- 
bulent flow dynamics. Viscosity, similarly as resistivity 
in a conductor, leads to entropy production by degrading 
inhomogencitics in the velocity field. While ideal fluids 
with ?y = cannot exist, it is interesting to seek for perfect 
fluids which come as close as possible to this ideal. 

Viscosity has the units of hn where n is some den- 
sity. To quantify the magnitude of the viscosity, it is 
natural to compare rj/h to the density of thermal excita- 
tions, nth, which can be estimated by the entropy density, 
s r*> ksrith- Motivated by the nearly perfect fluid behav- 
ior seen in the RHIC experiments, Kovtun et al. have 
recently postulated a lower bound for the ratio of r\ and 
s for a wide class of systems [5J: 
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Equality was obtained for an infinitely strongly coupled 
conformal field theory by mapping it to weakly coupled 
gravity using the AdS-CFT correspondence. While ex- 



amples violating the bound Qwere found (see [7]), the 
existence of some lower bound with ksTj/hs of order unity 
for a given family of fluids is not unexpected. It is analo- 
gous to the Mott-Ioffe-Regel limit for the minimum con- 
ductivity of poor metals [HI H], and to the saturation 
of the relaxation rate at — ksT/h ■ 0(1) close to 
strongly coupled quantum critical points |10j . In all these 
cases an exhaustion of scattering channels and a satura- 
tion of the kinetic coefficients occurs, once the mean free 
path becomes comparable to the interparticle distance. 
It follows, that the ratio 77/s is a unique indicator for 
how strongly the excitations in a fluid interact. 

While effects due to electron-electron interactions only 
amount to very small additive corrections to the conduc- 
tivity a{oj,T) in the collisionless optical regime hw ^> 
ksT [31 HI HI] , collisions are crucial in the opposite regime 
Hui <C ksT [21 [TS]. There they establish local equilib- 
rium, the remaining low frequency dynamics being gov- 
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FIG. 1: (a) Velocity profile u and associated Stokes force 
density f s = n\7 2 u counteracting the current flow, (b) Inho- 
mogeneous current flow expected in a four-contact geometry 
with split source and drain contacts held at voltage ±V/2. In 
the absence of viscous and other nonlocal effects, the current 
would be proportional to the applied voltage V, independent 
of the distance L between the contacts. Viscous effects di- 
minish the current as L decreases. 
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erned by hydrodynamics, i.e., by the slow diffusion of 
the densities of globally conserved quantities. Clearly, 
a hydrodynamic description works best when collisions 
are frequent and the fluid is strongly correlated. Trans- 
port in such a regime reveals information about the na- 
ture of the excitations. Examples ranging from the two 
He-isotopes [13] E], cold atomic gases [T5HTD] and elec- 
trons in semiconductors and metals [III [18], extending 
all the way to the high energy regime of the quark gluon 
plasma |5J and matter in the early universe |19j . 

In graphene, for energies below a few electron volts, the 
electronic properties are governed by the Hamiltonian 
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with the Fermi velocity v ~ 10 8 cm/s [2]. p/ = —ih\7 ri 
is the momentum operator, I = 1, N labels the N = 4 
spin and valley (2 Fermi points) degrees of freedom, and 
(7 = (a Xl <7 y ) are the Pauli matrices acting in the space 
of the two sub-lattices of the honeycomb lattice struc- 
ture. Without the Coulomb interaction, Eq. ^ is the 
Hamiltonian of N species of free massless Dirac par- 
ticles [2D]. The strength of the Coulomb interaction 
is characterized by the effective fine structure constant 
a — jf— ~ 2.2/e, which is not small for realistic values 
of the substrate dielectric constant e. Key for an un- 
derstanding of clean, undoped graphene is the fact that 
it is a 'quantum critical' system with marginally irrel- 
evant Coulomb interactions which renormalize logarith- 
mically to zero [3] [H [H] [21] [221 [231 [24] . This quantum 
critical behavior in undoped graphene is responsible for 
the distinctly different behavior of the collision-free and 
collision-dominated frequency regimes [25] , 

Collision-dominated transport can most efficiently be 
addressed by solving the Boltzmann transport equation 



d_ 
dt 



1 de 



kA 



1 de 



kA 



h dx 
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for the quasiparticle distribution function f — fx (k, x, i), 
which depends on momentum k, position x, time t and 
band index A = ± (labelling upper and lower parts of the 
Dirac cones centered at the two Fermi points). Eq. |3| 
can be derived from a nonequilibrium quantum many 
body approach, yielding the collision integral J^ co \i in 
terms of the Coulomb interaction |3J. The moments of 
the distribution function yield the conservation laws for 
the charge density, 9tp+V-j = 0, the momentum density, 



[<9 t u+(u- V)u]+Vp- 
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f s = 0, (4) 



and the energy density e. j is the current density and 
u the velocity field of the fluid with enthalpy density 
w = e+p where p is the pressure. For undoped graphene 
the Gibbs-Duhem relation implies furthermore w = Ts. 



Eq. Q is the Navier-Stokes equation for graphene, de- 
rived under the assumption |u| <C vf- Compared to non- 
relativistic hydrodynamics there is an extra relativistic 
term cx d t p, but at low frequencies its effect is small. 
The Stokes force /? = 9jT^ = ij\7 2 Uj is determined by 
the leading dissipative contribution to the stress tensor: 
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■qXij + C<5. y V ■ u. 
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in an expansion in gradients of u. Here = dui/dxj + 
duj/dxi — SijV ■ u corresponds to a pure shear flow while 
the second term is a volume compression. r\ and Q are 
the shear and bulk viscosity, respectively. 

In what follows we include all contributions to 3 m \\ up 



(2) to order a , keeping in mind that a flows to zero as T 



0. Assuming a slowly varying, divergence free velocity 
field u(r), we determine the shear viscosity by computing 
the stress tensor in linear response. Close to equilibrium, 
an inhomogcncous flow field constitutes a driving term 
in Eq. (|3j of the form v k>A • V x / A = £\ 



bijXji with 
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Here, vt,A = Vk£kA/^ is the velocity of quasipar- 
ticles with energy e^x, P — l/fegT and Iij (k) = 

{jW~ ~ l^ij)- I n li ncar response, the distribution 
function can be parametrized as: 
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where / cq = f\(]s.)\g ij= Q. Linearizing the Boltzmann 
equation in the zero frequency limit, it can be cast into 
an operator formulation: \<p) — C \g) [2DJ [27]. The op- 
erator C is Hermitian with respect to the inner prod- 
uct (a\b) = (8^ 2 ) _1 E uA / d2fca ^( k ' A )^( k ' A )- Thc 
gij (k,A) parametrize the non-equilibrium distribution 
function and are obtained by inverting the operator C. 
Using them to express the stress tensor and comparing 
with Eq. ([5| one obtains the shear viscosity: 
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The dominant contribution to tj comes from the small- 
est eigenvalues of C restricted to functions given by 
Eq. ([6]). The inversion of thc collision operator can be 
a formidable problem and usually requires a numerical 
solution. Thc problem simplifies, however, in two di- 
mensions where the amplitude for collinear scattering 
processes (involving quasiparticles with identical veloc- 
ity vector) is logarithmically divergent. Screening effects 
of higher order in a and lifetime effects cut off this di- 
vergence in the infrared at transverse momenta of order 
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aT/v [21 To logarithmic accuracy in a, we can thus 
consider collinear scattering processes only. The corre- 
sponding restricted operator possesses three zero modes: 
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which reflect the conservation of charge n, chirality x 
(the total number of particles and holes) and energy E 
in collinear processes. This conservation is exact for the 
modes g^ n,E \ while it holds only to lowest order in a for 
g^ x \ being due to kinetic constraints on the two-body 
scattering of massless Dirac particles (see [25] for a re- 
lated discussion). Eq. Q shows that these modes corre- 
spond to distribution functions which, when restricted 
to quasiparticles with identical velocity v = ue. re- 
duce to equilibria with direction dependent parameters 
/i(e), (p(e), T(e) conjugate to the conserved quantities. 

If there were only collinear scattering processes, the 
shear viscosity would be infinite. However, the inclusion 
of other processes causes rj to be finite. Nevertheless, 
the dominance of collinear scattering allows us - in lead- 
ing logarithmic approximation - to invert the operator 
C within the Hilbert space spanned by the modes (pL 
At zero doping, a divergence free velocity field does not 

(n) 

excite the mode g^ , and the relevant subspace is only 
two-dimensional. This inversion is easily done and yields 
{(f> |C _1 1 <fi) = C v 2~ 3 / 2 a~ 2 . The remaining numerical 
coefficient C v stems from the evaluation of the matrix el- 
ements of the full scattering operator in the 2d subspace 
of gfj' x and is expected to be of order unity. We obtain 
C v ~ 0.449, consistent with this expectation. The a~ 2 
dependence follows from the fact that the collision op- 
erator C is of second order in a. The shear viscosity of 
graphene finally results as 



v = c, 



N{k B TY 
Ahv 2 a 2 



1 + 



1 



log a 



(10) 



which is the central result of this Letter. It can be ra- 
tionalized by using the Fermi liquid result [T3] 7?fl — 
nmv 2 T rc \ with n — + ^thermal — (ksT/hv) 2 , relaxation rate 
t^j 1 ~ k B T/ (ha 2 ) [3 j and typical energy mv — > ksT. 
Extending the analysis beyond the leading approxima- 
tion by regularizing the logarithmic divergence in the for- 
ward scattering and including more modes gij , we obtain 
corrections of relative size l/log(l/a). For a = 0.1 they 



increase the leading result (10 1 by only 20% 



We need to keep in mind that the quasiparticle are 
not free, their dispersion reflecting the renormalization 
of the velocity v — > v [1 + f log(A/fc)] , where A is an 
appropriate UV scale. We implement this by a renor- 
malization group approach combined with scaling laws 
for physical observables. The coupling constant evolves 



as a (T) ~ 4/ log ^ (with T A = f^) while the velocity 
grows logarithmically v (T) = va/a(T). However, the 
combination [av] (T) entering r\ does not change under 
renormalization [3D]. Thus, Eq.(10l is the correct low 



temperature result for the renormalized quasiparticles. 

On the other hand, the entropy density of noninteract- 
ing graphene, including renormalization effects, is 4 
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The above finally results in the sought ratio: 
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As T — > the ratio r\js grows, a behavior expected 
for a weakly interacting system. However, since a is only 
marginally irrelevant rj/s grows only logarithmically. In 
contrast, in the regime T -C [i of doped graphene with a 
finite carrier density n, one obtains the usual behavior of 
a degenerate Fermi liquid [T31[TB] with 77 ~ hn(fi/T) 2 and 
s ~ kgnT/n, in which case r)/s ~ (h/kg) (/i/T) 3 di- 
verges much more strongly at low T, see Fig. [2] Higher 
order corrections of the long range Coulomb interac- 
tion [24] mainly reduce the regime where a (T) decreases 
logarithmically (they effectively reduce Ta). This fur- 
ther decreases the ratio rj/s. Similarly, additional short 
range interactions g yield leading additive corrections 
ex ag (T/T\) to the collision operator. Their effect is 
small provided that they do not lead to an excitonic in- 
sulator and the low energy physics of massless Dirac par- 
ticles is preserved [3"T] . 

Note the small numerical prefactor in (12 1. As shown 
in Fig. [2j it keeps the ratio rj/s small in a large temper- 
ature regime, where it approaches the value of Eq. ([l]). 
As was shown recently, cold atoms with diverging scat- 
tering length are materials which also come close to the 
value ([I]) ]T5][TS]. Our result shows that, interestingly, 
graphene has an even smaller ratio rj/s, thus being an 
even "more perfect" liquid than those critical systems. 
We envision interesting experimental manifestations of 
viscous effects in graphene, especially in the conductance 
properties of very clean samples. Viscous drag should re- 
sult in a decrease of the conductance with the linear size 
of geometries such as in Fig. [T] where two spatially sepa- 
rated contacts take the role of the source and the drain, 
respectively. Applying a source-drain bias V to these 
split contacts induces an inhomogeneous current flow and 
corresponding Stokes forces that oppose it. The smaller 
the spatial scale L, the larger the viscous forces, and 
hence one expects an increasing resistance. The latter 
would be scale-invariant in the absence of viscous forces 
and other non-local effects on conductivity [29] . 

The unusually low viscosity in graphene suggests the 
interesting possibility of electronic turbulence in this ma- 
terial. For simplicity we consider fluid velocities small 
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compared to the Fermi velocity, and analyze the low fre- 
quency limit of the Navier-Stokes Eq. Q. Turbulence 
arises from the nonlinearities oc (u • V) u while dissipa- 
tion due to the Stokes forces suppress turbulent flow for 
large rj. The dimcnsionlcss number determining the rela- 
tive strength of these two effects is the Reynolds number: 



Rc = £/*£ x ML x !%E 
rj/H hv/L v 



(13) 



where we assumed a typical fluid velocity u typ and a 
characteristic length scale L for the velocity gradients. 
Thus, the ratio i]/s reveals itself as the key character- 
istic determining the Reynolds number, apart from geo- 
metrical factors, typical energies and velocities. To ob- 
serve turbulence one needs Re > 10 3 ■ • • 10 4 in 3d, and 
somewhat higher values in 2d. However, even for lower 
Re ~ 10 • • • 10 2 two dimensional flow in the presence of 
extended defects or nano-sized obstacles undergoes com- 
plex phase locking phenomena and chaotic flow |32j . Ap- 
plying strong bias fields to graphenc, fluid velocities of 
the order w typ ~ O.lf can be achieved [331 . while still 
avoiding the onset of non-Ohmic effects. In the collision- 
dominated regime of undoped graphene the fluid veloc- 
ity induced by a field E scales like u/v ~ eEHv / '(fcgT) 2 . 
Hence, the Reynolds number increases as 1/T with de- 
creasing T, enhancing the tendency towards turbulence. 
With flow velocities as estimated above we expect com- 
plex fluid dynamics as in Ref. |32j already on small length 
scales of the order of L ~ 1/im. This would consti- 
tute a striking manifestation of the quantum criticality 
of graphene and could be relevant for potential nano- 
clcctronics applications of this exciting material. 
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FIG. 2: Ratio n/s in graphene as a function of T. The UV 
cut-off was taken to be Ta = 8.34 • 10 4 K following g]. In 
the undoped, quantum critical system n/s is very small in a 
large temperature window where the coupling a(T) remains 
of order O(l). The value 1/47T obtained for some strongly 
coupled critical theories is shown for comparison. As shown 
in the inset, away from zero doping and quantum criticality 
(T < \fJ,\), the viscosity assumes the behavior of a degenerate 
Fermi liquid, 7]/s ~ {\fJ-\/T) 3 . 
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